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Abstract 

Wc provide a constructive global existence proof for the multivari- 
ate viscous Burgers equation system denned on the whole space or on 
a domain isomorphic to the n-torus and with time horizon up to in- 
finity and C°°- data (satisfying some growth conditions if the problem 
is posed on the whole space). The proof is by a time discretized semi- 
explicit pcrturbative expansion in transformed coordinates where the 
convergence is guaranteed by certain a priori estimates. The scheme 
is useful in order to define computation schemes for related equation 
systems of fluid dynamics. 

2000 Mathematics Subject Classification. 35K40, 35Q30. 

1 Introduction 

In order to understand the dynamics of fluids we need to understand equa- 
tions of the form 



duj _ xr^n d 2 Uj _ v^n dv± 
dt - V 2^j=l -Qxf L>3=1 dxj 



J 



(1) 



I u(0,.) 



where v is some strictly positive constant (i.e. v > 0) and 1 < i < n. Global 
existence for such equations can be obtained by standard methods, e.g. by 
an upper bound estimate for the solution with respect to the maximum 
norm (cf. the estimate (|10p below). In this paper we look for a constructive 
global solution scheme for u = (u±,--- ,u n ) on some domain [0,oo) x 17, 
where [0, oo) is some time interval and O C M n , and where some initial data 
h = (hi,-- - ,h n ) T with hi G C7°°(ft) are given. Special interest is either in 
the case O = R n or in the case of periodic boundary conditions, i.e. solutions 
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of the form U{ : [0, oo) x R n — )■ R, 1 < i < n which satisfy the condition 
U{(t, x) = Ui(t, x + 7) for 7 £ T with some lattice F, say T = Z n . The latter 
periodic functions define solutions on the quotient space [0, 00) x R n /Z n in 
a natural way, and the spatial part of the latter is isomorphic to the n-torus 
T". Alternatively, one may look at the restriction of periodic functions of 
the kind above to a reference domain = £Iq, where 

Q = { x e R n \0 < Xj < 1, 1 < j < n} , (2) 

and then identify according to 

Ui(t, x) = Ui(t, x + 7) if x, x + 7 G <9f2o and 7 G r, (3) 

where, dQo denotes the boundary of Oo- We shall speak loosely of the 
Cauchy problem on the n-torus = T n . Clearly, the data h = (hi, • • • , h n ) T 
are assumed to be periodic if Q = T n . Our method can be adapted to 
other domains C W 1 with additional boundary conditions of Dirichlet, 
von Neumann or mixed type. Furthermore, in prinicpal our method can be 
applied to fluids on compact Riemannian manifolds M, where an equation 
analogous to ([1]) is of the form 

— -£u + V u ii = 0. (4) 

Here, V denotes the covariant derivative, and L may be the Hodge Laplacian 
A = -(dd* + dd*) or the Bochner Laplacian -V*V with V : C°° (M, TM) -> 
C°° (M, T* <g> T), or some similar operator arising from the covariant deriva- 
tive. 

However, in order to avoid too much technicalities we mainly stick to the 
case = M n , and confine ourselves to cursory remarks regarding Cauchy 
problems on = T n , on manifolds or initial-boundary value problems. Our 
goal in this paper is to define a constructive scheme for the multivariate 
Burgers equation which also can be interpreted in a probabilistic way. In a 
subsequent paper we shall show that such schemes may be useful in order 
to define global computation schemes for the incompressible Navier-Stokes 
equation. 

Back to the multivariate Burgers equation, if we put a potential source 
term on the right side of ([1]) , then we have a generalization of Burgers one- 
dimensional model for fluids, i.e. 

- z/Au+ (u- V)u = -VF(t,x). (5) 

Remark 1.1. If the velocity field equals the gradient of some potential func- 
tion (j> (in analogy to a physical potential of a conservative force) for initial 
time to 5 then we have 

u(to, x) = — V<f)(to, x) (to some initial time), (6) 
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and equation ([5]) is equivalent to the evolution of <f> with 



^-vA0=±\4>\ 2 +F(t,x). (7) 
So-called Hopf-Cole transformation may be used to transform to 

da . 1 „ 

— - uAa = —Fa. 8 
dt 2v y ' 

Global existence and effcient computation schemes for such scalar equations 
are well known and for some conditions on F explicit local expansions of 
the fundamental solution lead to schmes which are of particular efficiency 
(cf. [IQ]). We know that data u(0, .) satisfy ([6]) for some potential function 
(f) if and only if the curl of u is zero. Incompressibility is then given for 
free since V • (V x u) = div (V x u) = for smooth data. However, for 
studying fluids we need to analyze the situation for data more general than 
just potential data. We deal with this more general situation in this paper 
where no Hopf-Cole transformation is available. 

Remark 1.2. Note that global existence for the system ([1]) can be obtained 
from estimates of the solution in the supremum norm in terms of the supre- 
mum norm of the initial data, i.e. of the form 

maxsup \uj(t, x)| < maxsup \fj(x)\, (9) 
3 xen 3 xen 

which may be obtained form estimates of the form 
^\\u(t,.)\\ H s <\\u(t,.)\\ Hs+ iJ2 E \\D a u l D^u :i \\ L 2-2\\Vuf H s (10) 

i,3 \a\ + \/3\<s 

for some postive sGR. 

Systems of form ([TJ) are building blocks of well-known models in fluid 
mechanics and finance. Let us consider the Navier-Stokes equations 

— -i/Av+(v- V)v = -Vp + f t>0, xen, (11) 
at 

together with the incompressibility assumption 

V-v = 0, t>0, xel], (12) 

and together with the initial conditions 

v(0,x) = h(x), x G n, (13) 

where h(x) = (/ii(x),--- ,h n (x)) is a given function with components in 
C°° (f2). This equation is to be solved for an unknown velocity vector 
v(t,x) = (vi(t, x), ■ ■ ■ ,v n (t,x)) with (t,x) £ [0, oo) x fi, and v(i, x) £ 
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[C°°([0,oo) x n)] n , and with a scalar pressure p{t,x) G C°°([0,oo) x SY 1 ), 
where v = 2 > is a the strictly positive viscosity constant. Furthermore, 

f(M) = ,fn(t,x)) 

is a given external force where /, G C 00 ([0, oo) x $7) for each 1 < i < n. 

Remark 1.3. Let us mention here, that in case of a compact Riemannian 
manifold the Navier-Stokes equation takes the form 

- £u + V u u = -grad p, div u = 0, u(0, .) = h(.), (14) 

where we use the notation 'grad' instead of the 'nabla'-operator notation 
in this context in order to avoid confusion with the symbol of covariant 
derivative. On ker div one has 

£u = Au + 2Ricu, (15) 

where Ric denotes the Ricci tensor. For manifolds with constant curvature, 
i.e. Ric = c • id, we get L = A + 2c • id, and in this case we can eliminate 
the pressure and write (|14p in the form 

^-£u + PV u u = 0, u(0,.) = h(.), (16) 

where h(.) is still assumed to have divergence zero. Here, P is the Leray 
projection, i.e. the orthogonal projection of L 2 (M,TM) onto the kernel of 
the divergence operator. 

If = M. n , then a natural growth condition for the data is that for all 

x G R n 

(iTR 

holds for any multiindices a and K with some constants C a x- Furthermore, 
in this case it is assumed that 



|d a h(x)| < T, ■ u tf ( 17 ) 



\d?d a i(t,x)\ < {l ^ K +t)K for all (t,x) G [0,oo) x R n , (18) 

and for any multiindices a and and integers K with nonnegative integers m, 
and with some constants C am K- 

Remark 1.4. In the following we assume that f = 0, although the method 
considered is not restricted to this case. Note that the condition (fT7|) ensures 
that 

h G [H s (R n )] n , (19) 
where for each s£R H s denotes a Sobolev space, i.e. 

H s = H s (R n ) = {/€§' (M n ) |/ is a function and ||/||^ < oo} , (20) 
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along with 

11/112= / |/(2/)| 2 (l + M 2 r^<oo, (21) 

and where / denotes the Fourier transform of /, and §' (M n ) denotes a 
space of tempered distributions. We shall use these function spaces for the 
convergence of the semi-explicit perturbative expansion on the whole space 
[0,oo) x W 1 . In case p = 2 the definition is equivalent to the more classical 
definition of spaces H k ' p for nonnegative integers k and real p 6 [1, 00), i.e. 

H k,p ( R n) . = | u G L p ( R n) \ D a u £ L p ( R n) for < fc | _ ^ 

More generally, for s£R one defines 

H s,p ( R n) . = K -S]j> ( R n) ; where A s u = gr-1 ^ + y 2) S /2^ > (33) 

and 3~ _1 denotes the inverse of the Fourier transform operator. 

Next for an arbitrary step size T > in the original time coordinates we 
shall find a time discretization {7}, / = 0, • • • , N} with T\ = Ylk=i Pk t 00 
where we define a local scheme recursively for each time interval 
in such a way that the growth of the solution u is well-controlled in order to 
allow for a global scheme. We shall see that for our global scheme the step 
size is of the form 

Pi ~ j (24) 

for a fixed constant C such that the scheme is indeed global with respect to 
time. We noticed above that a growth condition is available from a priori 
estimates, but we shall get an estimate independently. One advantage of 
our approach is that our consideration can be generalized in order to de- 
fine global schemes for Navier-Stokes type equations. A second advantage 
is that we get a constructive solution scheme which leads to computation 
schemes. In order to prove local convergence with respect to time we shall 
consider a simple time transformation at each time step, i.e., we shall con- 
sider transformations of the form t = pir where t £ [7)_i,7]] , and where 
r € [(I— 1)T, IT] for natural numbers I > 1. On each such interval we denote 
kth approximation of the solution by 

(T,x)^u^ = «' fe 'V- - ,< M ) (25) 

where the index k in u p ' k ' 1 refers to iterations at a given time step and I — 1 
is the number of time steps accomplished. We start with the first time step 
solving (pQ) on [0, p\T] x W 1 with respect to t coordinates and on [0, T] with 
respect to r coordinates, where t = p\T (the parameter p\ will be determined 
below). Let the function (r, x) — > u p,1 '° = (u^' 1,0 , ■ ■ ■ ,Un' 1,0 ) be the solution 
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of the equation system 



Pi 



En d 
7=1- 



du 



p,l,0 



'3 dx j 



1 < i < n. 



(26) 



nP' 1 ' ^, .) = h. 



Since the hi, 1 < i < n are known, at this stage we have no coupling, and 
(|26p are essentially n scalar equations. Having defined u p ' k '° = (u1 ,k '°, ■ ■ ■ ,u n 
with u p ' k '°(0, x) = h for some k > 1 we define u /9,fc+1 '° recursively to be the 
solution of (corresponding to the domain [0, p\T] xl" in original coordinates 
t and to [0, p\T\ x W 1 in transformed coordinates r.) 



du 



P,fc+1,0 



ih 



pi I D,-=i - 



2„,p,fc+l,0 



En 
7=1 S 



p,k,0du 



P,k+1,0 



dxi 



1 < i < n, 



u p,fc+l,0( 0> .) = h. 

(27) 

Then subtracting the equation (|27p for u p,fe '° from the equation (|27p for 
u p,fc+i,o we g et 



dSu 



P,k+1,0 



dr 



Pi ELl 



p,k+l,0 



En 
.7 = 1 



p,fc,0 



p,fe+l,0 



piE,«' ^ 



(28) 



where <5tt: 



p,fc+l,0 



p,fe+l,0 

w7 



tif ,fc ' for fe > 1, and 



,1,0 



u 



p,l,0 



/ij. Now 



look at the family of equations in (|28p . They are linear equations with 
first order coefficients u p -' ,0 available externally from the previous step re- 
spectively. Again these are essentially n scalar equations. The coupling is 
'externalized' to the coefficients u? ,k '° and the source terms carrying infor- 
mation from the previous step. Next p\ > can be chosen such that on the 
domain [0, piT) x 1" 



(29) 



is a strongly convergent series and defines a function u p '° in a suitable clas- 
sical space. Having defined the scheme for time steps 0, • • • , I — 1 we define 
the scheme with time step index I (i.e. the scheme for the solution on the do- 
main [T{_i, T{\ in original time-coordinates t) in an analogous way (using the 



semi-group property). Let the function (r, x) — > u 



p,i,l 



« 



p,i,i 
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be the solution of the equation system 



dw 



p,i,i 



pi £i=i-a£3- -£,•=!«: 



dxf 



1 < i < n. 



(30) 



u^((Z-l)T,.) = u^- 1 (G-l)T,.), 



and let (t,x) — >• u p ' fc+1 '' = , • • • ,u„ fc+1 ' Z ) be the solution of the 

equation system 



Ou 



p,k+l,l 



Pi I 2^=1 —fa 2^=1 



n p,k,ldi£ 
UY — - 



P,fe+i,i 



1 < j < n, 



u p,k+i,i((i _ .) = u P,i-i((2 _ i)r, .). 

(31) 

Remark 1.5. Note that in the original coordinates the initial condition at 
the Ith time step is 

u 1 ' , ((T,_i,.) = u"' 1 ''- 1 ((Z-l)r j .). 
Then subtracting the equation (|3ip for u p ' fc,i from the equation (|3ip for 



95« 



9t 



9 8u- 

pi E,-,i — &j E,-,i »' 



n o.k.ldSu 1 ^ 



p,fc + l,i 



5u p 



k+l,l 



0. 



(32) 



This leads to a functional series 



p,l p,i,l 
uy = uy 



p,k+i,i 



\<i<n 



(33) 



which is defined on [(/ — 1)T, IT] x lR n in r-coordinates corresponding to 
[Tj_i,Tj] x W 1 in original t-coordinates. If we can choose pi, I > 1 such that 



(34) 



then we have a global scheme. We shall prove in this paper that such a 
choice is possible such that the scheme is indeed global. Moreover, we shall 
see that the convergence of (f2l?|) and (|3"3"j) is such that differentiation can be 
done term by term, and convergence of the differentiated functional series is 
uniform and absolute in a pointwise sense. Plugging in the series 
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into (H]) and rearranging the summands according to the laws one checks 
easily that the limit in (|29l [33]) leads indeed to a solution of the multivariate 
Burgers equation. 

The construction outlined has an interesting consequence: the solution 
of the multivariate Burgers equation has a representation in terms of a series 
of fundamental solutions of (scalar !) linear equations. This may be use- 
ful for numerical purposes since there are efficient schemes for fundamental 
solutions of scalar linear parabolic equations in terms of WKB-expansions 
or related analytic expansions. Note also that this leads to probabilistic 
schemes where elaborated weighted Monte-Carlo schemes may be used (cf. 
[7] and [S] and references therein). This would be attractive especially in 
cases of complicated boundary conditions, or if coupling with other equa- 
tions (e.g. Maxwell equations) lead to complcated models where traditional 
techniques are of very limited success so far. However, the following semi- 
explicit formulas may also be used in the context of sparse grids or adaptive 
sparse grids etc. (regularity for the solution makes this option attractive). 
Let us spell out the scheme in terms of fundamental solutions (densities). 
For each positive integer k > 1 and given 1 < k < n and time step I 
let us denote the fundamental solution of the linear equation 

dwi _ / v^n d 2 Wj , v^n .,p,k,l du>j \ (nr\ 
dr - Pi \L>j=l ~dxf + L>3=1 U j fej J [ ' 

by T^ ( we can construct this solution by the classical Levy expansions 
and find very effcient higher order approximations in terms of local ana- 
lytic expansions; the upperscript f2 of r£' n,J indicates that we apply such 
approximations on bounded domains obtained either by transformations or 
by cutoff). Then given u p ' k ' 1 we have the representation 

6u?' k+1 ' (t,x) = 

(36) 

Jo In P £j Suf k ' l (a, y)Jffc—(a, y)T p k ' n ' l (r, x, a, y)dady. 
The series starts at I = with the fundamental solution T^' '° of 

^ = p(pUi% + £?=i^t|)> (37) 

and with 

Ju^'V^) = Uj' 1 ~ hj = J hj(y)T p ^ a ' (r,x;0,y)dy - hj(x), (38) 
and proceeds analogously for I > 0. We can make the representation fj36[) 
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at each time step I more explicit, i.e. 
Suf +1,1 (t,x) = 



lo In'" lo" In ( E£=i Pi 5u h' l ^£-) i^Vi) 

-rrk I V^n p,mldu f ' ,m ' 1 \ , v 



Il^r^'V™. Vm, cr m _i, y m _i)r^' ' (r, a kl y k )daidyi ■ ■ ■ da k dy k . 

(39) 

Hence we have the formal representation for the solution u = (u\, • • • , u n ) T 
of the multivariate Burgers equation with zero source term 



< 



• 1 '°+EEo (Er=oK ,fc+1, fa*) 

<' 1>0 + E£o /o T In" lo" In fc 1= i ) Vi) 



Eoo -rrk I \-^n p,rn,l du?' m ' 1 \ , v 

fc =0 n m=2 (Ei m =lA>i«£, J (^m,ym)x 

n m= X 2 r m ''(o- m ,y m ,crm-i,2/m-l)r^ ' (r,x,a k ,y k )daidyi ■ --da k dy k , 

(40) 

where for each I the whole related series for k has to be summed up, and 
where 

u^ 1,0 (t,x) = I a hj(y)T{' ' (T,x;0,y)dy, 1 < i < n 

£- k uf' (r,x) = f a h j (y)^Tf l ' (T 1 x;0,y)dy, l<i,k<n (41) 

Suj' '°(r, x) = Uj' ,0 (r, x) — hj(x), 1 < j < n. 

Analogously for time steps I > 0. The latter recursive scheme leads to 
semi-explicit recursive formulas for the multivariate Burgers equation. Note 
that the fundamental solutions of the parabolic recursively linear equations 
involved have an explicit representation. It leads also to new higher order 
numerical schemes. 

In the next section we show how a series (pi) can be chosen such that we 
have global convergence of the scheme above. 
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2 Convergence of the global scheme 



We may assume that the time step size in transformed coordinates equals 1, 
i.e. w.l.o.g. we assume T = 1 (the scheme of section 1 is a global scheme if 
it converges to the solution, where Y2i Pi t °°)- This leads to a time scheme 
with time-step size of the Ith time step proportional to j . 

Remark 2.1. Justification of the definite article in 'the solution', i.e. unique- 
ness, is provided below in Theorem 3.1. 

We construct the solution u(t,x) = u p (t,x) of the multivariate Burgers 
equation in terms of a time discretized scheme of functional series where for 
all Z > 1 and for (t, x) G [Tj_i, T{\ x W 1 resp. (r, x) G [(/ - 1), I] x R n we have 

<(r, x) := u p /(t, x) = uf' l (r, x) + 5uf +l \r, x), 1 < i < n. (42) 

k 

We show that there is a sequence of numbers (/>&) and a related time dis- 
cretization {T[\l = 1, • ■ • , N} where T\r f oo as N f oo such that u p = 
(u^, • • • ,Un) T defined via ([12]) satisfies the multivariate Burgers equation. 

Remark 2.2. It is interesting that the condition of decay for the initial data 
in ()17p implies that we can transform the Cauchy problem to a compact 
domain considering the transformation yi = arctan(xj), 1 < i < n and the 
related equation for the function u, where 

Ui(t,y) = Ui(t,x) resp. u^(r,y) = u?(t,x). (43) 

Then we have 



(44) 



dxj dyj 1 + tan 2 (?/j) ' 
and 

d 2 Ui d 2 Ui 1 diii 2tan(y^) 

d Xj dx k = 8y] (l + tan 2 ( yj )) 2 " jk Wj (1 + tan 2 ( % )) (1 + tan 2 (y fc )) ' 

_(45) 

where <5jj denotes the Kronecker delta. The Cauchy problem for u = 
(u±, ■ ■ ■ u n ) T related to the multivariate Burgers equation live on the domain 
[0, T] x (— 7r/2, ir/2) n . It is interesting that, for example, the relation 

|2 ^3 d 2 Uj 

Cjdx k 

enforces 



(l + |x| 2 ) 3 ' |0as|x|too (46) 

v ' OX^OXi 



02- 



|0 as |x| tvr/ 2 , ( 4 7) 



dxjdxk 

and our analysis below shows that ()46p is indeed the case if the initial data 
satisfying (|17|) . This makes Dirichlet data natural for the related initial- 
boundary value problem for u. We shall use this in the second part of this 
paper when we deal with the Navier Stokes equation. 
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In order to show the global convergence of the functional series u% we use 
classical results for scalar equations. First we shall consider convergence with 
respect to the supremum norm | . | q , where for bounded mesurable functions 
we define 

\g\o-= sup \g(t,x)\. (48) 

(t,x)e[0,T]xK™ 

As we shall see the growth of the solution u of our scheme is controlled 
by the maximum principle. This growth is linear with respect to the time 
horizon and controls the time step size of the scheme (which has reciprocal 
value essentially). In order to ensure the convergence of the substeps k of 
the Ith time step (which lead to a local solution of the multivariate Burgers 
equation with respect to time) we observe that we can control the first order 
derivatives of the correction terms 5u? ,k ' 1 uniformly. The time step size is 
chosen in such a way that we get a uniform and absolute estimates with 
respect to the norm |.|i2- At each time step we also establish an uniform 
and abolute convergent geometric series bound for the correction terms series 
(^u^' k ' l ^j for each 1 < i < n with respect to the norm |.|i,2- This implies 

that derivatives can be taken term by term and that the functional series 
(|42p is indeed a solution of the multivariate Burgers equation. So much for 
the outline. 

Next we need some elementary technical preparations (quite standard). 
Define the Euclidean distance in M n+1 between the points y\ = (ti, xx), y 2 = 
(t 2 ,x 2 ) by 

e(z 1 ,z 2 ) = \/\t\ - t 2 \ + \xi - x 2 \. (49) 

Although we do not consider convergence w.r.t. Holder norms we shall use 
them at one step of the argument. If w is a function in a domain D C IR n+1 
we denote 

r i \w(yi) - w(y 2 )\ 
[w}s/2,5,d = , sup . (50) 



Next define 



\Ms/l,S;D = \w\o,D + [w\s/2,S;D; (51) 

Furthermore we shall use some classical spaces 

n 

IHo,l;D := \w\o;D + ^ \w Xl \ Q - D , and (52) 



|^ll,2;J) — \w\0;D + y) \Wxj\o;D + \wt\o-,D + 'Y] \w XlXj \o-D- (53) 

Note that the latter norms do not define Banach spaces. However, we 
shall use the fact that a functional series which is uniformly and absolutely 
bounded with uniformly and absolutely bounded derivatives can be differ- 
entiated term by term. For this matter (|53|) is useful. In the following 
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for T G (0,oo) let D = Q T = [0,T] x IET. For vector-valued functions 
w = (wi, • • • , w n ) T we define 

|w|iq£) = max |i0i|i,2.D (54) 

l<i<n 

etc. Note that we can rewrite the Cauchy problem (JT]) in terms of coordinates 
r(t) = for r 6 [Tj-i, T{\ as in ([58|) where we define 

p(r) if r€ [Tn.Tf] (55) 

for all I > 0. Well, on key idea of the global scheme of this article is that 
the time step size decreases with order j at the Ith time step where the 
estimate of the solution is bounded by C\ = Co + 1 with Co a bound for the 
initial data (with respect to some strong Banach norm). This implies that 
the first order coefficients at the iteration substeps of each time step are 
uniformly bounded by some constant C* which is independent of the time 
step number I. 

Remark 2.3. Note that a global (non-constructive) estimate (as ([U]) based on 
()10|) above in remark 1.2) may lead us to a global scheme with an uniform 
time grid. However, we do not consider this for two reasons. First we want 
an explicit control of the time step in terms of the data and this is not 
provided by the a priori estimate (the size of the time steps may be uniform 
but so small that they are not useful in practice). Second, the present 
scheme may be extended to Navier-Stokes equations and there are no such a 
priori estimates available for the Navier-Stokes equation in the crucial case 
of dimension 3. 

We have 



Theorem 2.4. Assume that the initial data h satisfy the condition (11). 
Let 

C := lu^lo := max |<' 1,0 | , (56) 

j£{l,-n} J 

and define C* as in lemma [2751 below. Define positive real numbers (pi) and 
(Ci) recursively via 

Pi = 4C*CVi = 4C*(C + (i-l)) ( 57 ) 

where C\ = Q_i + 1 for I > 1. For this sequence (p^) the functional series 
scheme above converges to the global classical solution u p = (u±, ■ ■ ■ ,Un) T €E 
[^([0,oo)x(l)]" of 



Qui _ t \( d 2 u? P^ u i\ 



(58) 



h. 



where (t,x) — > u(t,x) = u p (t,x) solves the original Cauchy problem FT]) on 
[0,oo) x M n . The Ci control the growth of the solution with respect to time. 
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Proof. Recall that T = 1. First we start with the scheme of the functional 
series (|29p and prove its strong convergence with respect to the norm (|54p 
on the domain [0,pi] x M n = [0, T±] x R" in original time coordinates (cor- 
responding to the domain [0, 1] x M. n in transformed time coordinates (i.e., 
r). Note that p(r) = p\ on this domain. 

Since the n equations (for 1 < i < n) in (j26|) are n identical scalar 
equations and the initial data and first order coefficient functions hi are 
smooth and have bounded derivatives, classical results tell us that a unique 



smooth solution in C°°, i.e., (r, x) — > u 



(0,1,0 



« 



p,l,0 



,P,1,0> 



of the linear 



parabolic equation (|26|) exists which is bounded with bounded derivatives 
where all the iif' 1,0 are equal entries. Hence, by the maximum principle (cf. 
part (i) of lemma 1231 below) for Cq := maxjg^... „} supr T)X \ e Q 1 |/tj| > we 
have 

max K' 1,0 | < C . (59) 

i6{l,— ,n} 

Moreover, according to part (ii) and part (iii) of lemma 12.51 there exist a 
constant C* > such that 



sup 

(-r,z)eQi 



^ P,1,0/ \ 

(r ' x) 



< C*Co 



sup 

(t,x)gQi 



- 1,0/ \ 

^/ (r,x) 



at 



< C*C 



(note the original t- variable here), and for all 

o 2 



sup 

(t,x)£Qi 



p,l,0 



(r,x) 



< C*Q 



o- 



Define C* = (2 + n + n 2 ) C*C*C* for convenience. Choose 

1 



(60) 
(61) 

(62) 



Pi - 



4C Q' 



and assume inductively that 



max 

ie{i,- ,«} 



o,i •" 



max 



i,ie{i,-,n} 



9 





p,fc,0 



p,l,0i 



< max —r 

o i<j<n 2 fe 



9a; i 



5u p,i,o 



(63) 



(64) 



(65) 



First we verfiy that the relations (|64p and (j65[) are inherited at step + 1. 
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We have 



|<5u" fc+1 '°|o < Pl|^ fc '°| (C*C + En 



jm=2 



dxj 



dxj 



(66) 



< \6uP> k >°\o( Pl C*Co + pi 
<\5u^°\ (l + Pl C*C ) 

< h\Su^\ < 5 fft|<5u p ' 1, ° 



In order to close the induction we need estimates of the first order deriva- 
tives because these are involved in the source terms on the right side of the 
recursively defined equations. We shall use the fundamental solution. In 
the following we understand the term 'fundamental representation' in the 
following sense: if (t,x) — > f(t,x) is a function on a domain [0,T] x R n and 
(t, x; s, y) — > p(t, x; s, y) is the fundamental solution of a parabolic equation 
then we call the function 



(t,x) -> (f*p)(t,x) := 



f(s,y)p(t,x;s,y)dsdy 



(67) 



a fundamental representation (we have classical convolution if p is a normal 
density and is dependent on x — y only). The solution 5u?' 2 ' (for 1 < i < n) 
and its first order spatial derivatives can be represented by fundamental rep- 
resentations involving the source term of (|28p and the fundamental solution 
of (|26p . or the respective derivatives of this fundamental solution. Hence, 
according to lemma [231 and remark T2.6I below we get 



IK 



2,0 I 



max i,me{l,-" ,n} 

< pi maxiefi,... >n} YTj=i 



dx m 0U i 



(68) 



< piConC^Su^'X < Pl C*\5uP^°\ 



where r^' 1 denotes the fundamental solution of (|26|) (which is identical for 
each i). Here we use lemma [231 again. Note that we have |<5u p ' 1,0 |o on the 
right side of the latter estimate and not the derivative. Next we assume 
inductively that for m < k we have 



\6urX < (piC* n C Q 



m—1 



\5u 



p,l,0i 



(69) 
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Then we get (a rough estimate is sufficient) 



a<5 U p,i,o 



dx-j 



PiC*Cq 



<\^ P ' k '°\o{piC*C + PlyT ^ U5U - o 



yk 



m=2 



QS U P< m '° 



(70) 



Note that 

ip p,l,0i I p,l,0 , I 

\5u^ |o = - hj\ 



PL 

sup 

s,x 



9 1,0/ \ 
m «j (s,x) 



ds 



Q < C*C 0P1 < \. 



Hence with the choice 



(71) 

the functional series (j29|) is absolutely and uniformly convergent with re- 
spect to the norm |.|o,i- Note that the relation (17) concerning the decay 
of the initial data at spatial infinity is preserved for the all functions u^' k, °. 
Applying lemma 12.51 once more by a similar argument we get the absolute 
and uniform convergence of the first time and the spatial derivatives of sec- 
ond order. Hence, the first order derivative with respect to time of the series 
([29)) can be computed componentwise, i.e. 



9 n,o d p i q ^ d p fc+i,o 1 , ■ , 



(73) 



and the same is true for the first and second order derivatives with respect 
to the spatial variables, i.e. for all 1 < I < n we have 



8 



-u 



p,0 



d p ,i,o _,_ 9 Su p,k+ifi 



-u 



OX\ OX\ 

and for all 1 < I, m < n we have 



dxi 



1 < i < n, 



(74) 



2 



dxidx m 



-It' 



p,0 



o 2 



dxidx m 



p,l,0 . 



2 



dxidx m 



-Su£ 



p,fc+l,0 



1 < i < n. 



(75) 



Using these relations we observe that the function u^' 1 defined in terms of 
the series in (|29|) satisfies the multivariate Burgers equation in a classical 
pointwise sense. Hence, we have constructed a solution on [0, T\] x W 1 in 
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original time coordinates (reps, on [0, T\\ x W 1 in r- coordinates). On this 
domain we have 



K lo = \u 



f>° + Y,6uf +1 X<Co + l, (76) 



fc>i 

and 



+ ^8uf +1 X 2 <C* n (C + l). (77) 



,P,1|, „ _ |„P,1,0 

fe>i 

Using this estimate we get a lower bound for the next time step. We define 

92 = AC*C[ = 4Q(Co + l)' (78) 

As it turns out this choice controls the first order coefficient in the equa- 
tions of the second time step such that C* is again a uniform estimate for 
the convolution with the fundamental solutions T^ k . Let us observe this 
inductively and more closely. Assume inductively that we have defined Q_i 
(up to some I > 1), and 

1 _ 1 

91 ~ 4C*(C,_! + 1) " 4C*(C + /-1)' ( 9) 

and that for I — 1 > 1 the function u p ' i_1 has been constructed on the domain 
[(1 — 2), (I — 1)] x R n (corresponding to the domain [l]_2,I}_i] xR" in original 
time coordinates). We want to construct the solution \i P)l on [(I — 1), /] x W L 
(corresponding to the domain [T{_i, ?}] x M n in original coordinates) via the 
series 

< l = < xl + E 6 < k+hl ( 8 °) 

k>\ 

for 1 < i < n. 

First, we consider the equation (|30p (recall that T = 1). Similar as in the 
first time step we observe that the n equations (for 1 < i < n) in (|3U|) are 
n identical scalar equations and the initial data and first order coefficient 
functions u^ 1 ' 1 are smooth and bounded. Classical results tell us that a 

unique smooth solution (i.e. in C°°) (r, x) — > u p,1,/ = (u^ 1,1 , ■ ■ ■ , v^ 1 ' 1 ) exists 
which is bounded with bounded first order time derivative and bounded 
spatial derivatives up to second order, and which has equal entries. From 
the previous step we have an estimate of the first order coefficients in (|30p . 
i.e. 

pi max M^ l -\ 2 <piClCi<\. (81) 

The choice of pi implies that we can choose the same constant C* as before 
(applying lemma [231 as described in reamrk T2.6I below) in order to get 

max \ul' l \ 1>2 < C*Q. (82) 
ie{i,-,n} 
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Applying lemma [231 we get 



\5u<>> 2 ' l \ h2 < Pl max i6{lr . , n} 5u^^- * r< 



UK' 1 " .ri 



< W C*Q|<5u^| 



1,2 (83) 



where 1^' denotes the fundamental solution of f|30j) (which is identical for 
each i) and the constant C* can be chosen the same as in the previous time 
steps (recall uniform upper bound for the first order coefficients implied by 
choice of time step size pi). For the choice of C* we use again lemma 12.51 as 
described in remark l2~lfl We may assume w.l.o.g. that C\ > 1. Now assume 
inductively (we look at inductive substeps now) that for 2 < m < k we have 



max \Suf' l \ 1>2 < (ptQC*)™- 1 \5u^\. (84) 
ie{l,-,n} 



Then we get 



\Su p ' k+1 \ 2 < P i\Su p ' k ' l \ (c*Ci + £ 



d5u.P' m ' 1 



sk 

<\5^ k \(\ + Pl C* n C ) (85) 

Here we observe inductively that C* is a uniform constant for estimates of 
the convolution with T^f and its first order derivatives with respect to time 
and first and second order derivatives with respect to space. Hence we have 
convergence of the functional series (|80p to the norm |.|i.2- Again, uniform 
and absolute convergence of derivatives (first order with respect to time and 
up to second order with respect to space) with polynomial decay at spatial 
infinity leads to the conclusion that the functional series (|80p can be derived 
term by term. Hence, the first order derivative with respect to time of the 
series (|80p can be computed componentwise, i.e. 



d_ P ,i _ d_ P ,i,i \ ^ _9 

dt u% ~ dt Ul + ^dt 



i,i . sr^ " x p,k+i,i , , . , 



and the same is true for the first and second order derivatives with respect 
to the spatial variables, i.e. for all 1 < / < n we have 

9 P ,i d p i i & x p,k+i,l -, ^ ■ fQT\ 

k 

and for all 1 < /, m < n we have 

< l = + T, J^rH' k+hl , 1 < i < (88) 



dxidx m % dxidx m 1 ^ dx\dx 

k 
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Using these relations we observe that the function u p ' 1 defined in terms of 
the series in (|80p satisfies the multivariate Burgers equation in a classical 
poiintwise sense. Hence, we have constructed a solution on [0,7]] x W 1 . On 
this domain we have 

/c>i fc>i 11 n 1 

Using this estimate we get a lower bound for the next time step. We define 

Pl+1 = 4Q(Q + 1) (90) 
Hence, the sequence (pi) satisfies 

pl+1 ~ iC*(C + l) (91) 

such that ^4=1 Pi t 00 as N f oo. This implies that the scheme is global in 
time. 

□ 

In the following we state the lemma used above in order to get the 
uniform constant C* which depends only on the dimension n and the Holder 
norm of the first order coefficients of the subproblems. This Holder norm is 
uniformly bounded due to the choice of the 'time step size' p\ (especially, it 
is independent of the substep k and independent of the time step I). The 
independence of / is in part due to the fact that (with the choice of pi and 
Ci in the theorem) we have piCi < C for some constant C independent of 
I. Recall that C\ is the |.|o-bound of the Cauchy problem of the Ith time 
step. The independence of the substep number k is due to the fact that 
there is a uniform bound of the first order coefficients for each subproblem 
of the iteration. The constant C* is determined according to the following 
lemma (in the remark below we explain in further detail how it is used for 
convenience of the reader). We need the following facts: at each time step 
I we solve the multivariate Burgers equation on the domain [T}_i,T/] x K" 
(original domain) . The final data of the time step I are the initial data of 
the time step / + 1 and figure as the first order coefficient functions of the 
first substep problem at time step I + 1. For this reason we have to ensure 
that the initial data of the next time step are Holder continuous (in order to 
get a classical solution for the linear problem of the first substep of time step 
/ + 1). Note that we encounter linear parabolic equations in substeps. In 
the A;th sunstep of time step / + 1 the corection terms 5u p ' k,l+1 , k > 1 have 
a fundamental representation with source terms. In order to estimate these 
representations we need the source term to be bounded continuous. Since 
the source terms involve the first derivative of the solution of the last time 
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step I, we need bounded continuity of the first derivative to be uniformly 
controlled (independent of the time step) . The following lemma and remark 
summarize the results needed and how they are applied in our situation. 

Lemma 2.5. Let T > be some horizon. Consider the equation (recall that 
< v < oo is some constant) 

= uAu + Y / b l (t,x)— + g (92) 
t=i 1 

along with the initial condition 

u(0,x) = f(x) (93) 
on the domain D := [0, T] x R n . Then we have: 

(i) Assume that u is bounded and continuous, that the time derivative Ut 
exists for any t G (0, T), and that for any t G [0, T] the spatial deriva- 
tives u Xi and u XiX . exist and are continuous. Furthermore assume that 
the coefficients bi are bounded on D for 1 < i < n, and that u satisfies 
the equations $9ty) and [93) . Then we have 

Ho <T\g\ + \f\ . (94) 

(ii) Assume that T = 1. Let g = 0, and assume that |/|i < Cj, and 
\b\o — Cft? Ho,i < Cl holds for finite constants Cj,C^,C^. Then the 
Cauchy problem \92) . \93) has a unique classical solution u G C^' 2 , 
and for some constant C* > ( which depends on dimension n and on 
v only) we have 

Mo,i < (i + c?c?)|/|i 

(Hi) Let g = 0, and assume that |/|2 < Cf, and | b| o,i < Hi,2 < C 2 
holds for finite constants Cf,C},C?. Then the Cauchy problem \92) . 
h93\) has a unique classical solution u G C^' 2 , and /or some constant 
C| > (which depends on dimension n and on v only) we have 

Hl,2< (1 + ^C 2 *)I/|2 

(raj Assume that T = 1. Let / = 0, and assume that \g\ a /2. a — C g , 
and |6|i,2 < Cfe holds for finite constants C g ,Cb. Then the Cauchy 
problem \92) . \93\) has a unique classical solution u G C^' 2 , and /or 
some constant C** > (which depends on dimension n and on v, and 
on the Holder constant C g only) we have 

Ml,2 < 
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Proof. Assertion (i) follows from the maximum principle. Next we prove 
assertion (ii). Since v > is constant the b{ are Holder continuous with 
respect to the spatial variables uniformly in t the fundamental solution p 
of ([92]) exists. Since / S C| classical theory tells us that a unique classical 
solution u exists on [0, T] x M n , and (in case g = 0) u has the representation 

u{t,x) = [ f(y)p(t,x;0,y)dy. (95) 

JTSP 

Recall that the fundamental solution p in the Levy expansion form is given 
by 

p(t,x;s,y) := N(t,x;s,y) +/ / N(t,x;a,£)(f)(cr,£;s,y)dcrd€, (96) 



where 

N(t t x; S , „) = exp (-^f ) , (97) 

and <^> is a recursively defined function which is Holder continous in x, i.e., 

oo 

<f,(t, x- s,y) = J2 (LN) m (t, x; s, y), (98) 

m=l 

along with the recursion 

(LJV)i(*,a;;s,j/) = LN(t,x;s,y) 

(99) 

(LiV) m+ i := J n LN(t, x; a, £) {LN) m (a, £; a, y)^. 

Now, since the time derivative applied to iV cancels out with the Laplacian 
applied to N we have 

LN =iz h ^.- ( io °) 

j=l 3 

Since v is a constant the iV depends on x — y and t — s. Therefore we write 
N*(t — s,x — y) := N(t,x;s,y). Hence, for the spatial derivatives of first 
order we get 

J* R n f{y)p Xj (t, x; 0, y)dy = f Rn f{y)N*. (t, x - y; 0)dy 
+ Jr« Jo Jm™ ^ (*' x ~ & ^M "' °> V)f{y)dadidy 

(101) 

= - J R n f Xj (y)N*(t, x-y; 0)dy 

~ Jr» Jo J r ^* (*> a; - f ; cr)$(<7, f ; 0, y)/^ (y)dad^dy 
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by partial integration and change of order of integration. Here 

oo 

$(t, x; s,y) = J2 ( KN *)m(t, x; s, y), (102) 

m=l 

with 

(KN*) 1 (t,x;s,y)=KN(t,x;8,y) 

(103) 

(KN*) m+1 := J Q (LN*) m (t, x- a, £) (KN*), (a, £; a, y)d<7#. 
Here, 

n 

KiV(t, a; ; S ,y)=^6 i (i,y)iV*(i- S ,x-y). (104) 
i=i 

and (LN*) m is defined as above. The estimate of part (ii) of lemma I2~31 then 
follows from (jlOip and the standard estimates 

N*(t-8,x-y) < Q - 



(105) 

-^Ej s i x ~y) < ' 

where < /xo < 1, and 1/2 < /Ui < 1. Similarly, for the spatial derivatives 
of second order we get 

J R n fxjx m (y)N*(t, x-y; 0)dy 

+ /r» Jo Jr« ^ (*> * - 6 ^ofo f; 0, !/)/*, (y)dad^dy (106) 

+ / R „ Jo /k» ^ (*> x _ 6 CJ ) $ i( CJ ' f I °. l/)/*j«m (y)dcrdCdy. 
by partial integration and change of order of integration. Here 

oo 

x; s,y)=J2 (K*N*) m (t, x- s, y), (107) 

m=l 

(K*N*) 1 {t,x;s,y)= K^N(t,x;s,y) 

(K*N*) m+1 (t,x;s,y) := £ J u K^t,x; p, V ) (LN*) m (p,r],x;a,ti)x 

(K*N*) 1 {a,C;s,y)dad^dr ] dp, 

(108) 

where 



Kt(t,x;p,ri) = £ (^(/^)) 



P.z-tj), (109) 



21 



and 

$5; = 2$. (110) 

It is worthwile to consider the lowest order term of the latter Levy-type 
expansion. Iterated partial integration leads to the representation 

u(t, x) X[Xm (t,x) = J Rn f XjXm (y)N* (t,x-y; 0)dy 

+ / R „ £ /*« K (t, x-z;a)( (E, (a^&j (t, x)N(a, z - y)) f Xj (y)) + 
(2 Zj (b, (t, x)N(a, z-y))) f XjXm (j,)) dad£dy 
+ terms with weaker singularity integrands 

(111) 

Note that representations of the first order time derivative can be expressed 
in terms of spatial derivative up to second order. Hence the estimate in 
part (iii) of lemma 12.51 follows from the representation in equation (jl06D 
with lowest order term of form (jllip . Finally the estimate of part (iv) is 
standard (cf. [3]). 

□ 

Remark 2.6. At each time step the lemma [231 is applied with T = 1 in trans- 
formed r coordinates. Assume that / time steps have been performed leading 
to a solution u p ' 1 which is defined on the time interval [0, /] with respect to 
transformed coordinates r. Then part (i) of lemma l2.5l (i). i.e. the maximum 
priciple is applied in order to show that the supremum norm of the function 
bounded by the supremum of the initial data of the (Z + l)th time 
step, i.e. the quantity \u p ' l (l, .)\q. Then part (i) of lemma 12.51 with zero ini- 
tial data and recursively defined right side is applied recursively. The choice 
of the time step size pi ensures convergence of the functional series with the 
elements 5u p,k ' l+l for k > 2 with respect to the supremum norm. In order to 
construct a classical solution we need estimates of the first and second spatial 
derivatives and of the first time derivatives. In order to estimate |ii p,1, ' +1 |o,i 
part (ii) of lemma [231 is applied, and in order to estimate |w p,1 '' +1 |i,2 part 
(iii) of lemma 12.51 is applied. From the proof of part (ii) of lemma 12.51 we 
see that the quantity bounded by the the quantity \u p ' l (l, .)|i 

plus a series involving first order coefficients bi and multiple integrals of first 
order coefficients bi times the quantity \u p,l (l, .)\\. Since the first order co- 
efficients have a multiplier pi from the representation in the proof of lemma 
12.51 we observe that the choice of p\ of form iC \ c with some C* depend- 
ing only on v and the dimension n (which makes the first order coefficient 
of the Cauchy problem for u p ' 1,l+1 small) ensures that we have the bound 

K' 1,/+1 |o,i ^ \u p ' l (l, Oil U + C i-p lC u ) for some constant C > 0. Note that 
powers of p\ come from the representation of the higher order terms in Levy 
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expansion and are estimated here in form of a geometric series. Similarly, 
from the proof of part (hi) of lemma [231 we see that the quantity \u p > 1 ' |i 2 
is bounded by the the quantity \u p ' 1 (I , .)\2 plus a series involving first or- 
der coefficients bi and first order derivatives of bi and multiple integrals of 
first order coefficients bi and its first order derivatives times the quantity 
\u p > l (l,.)\ 2 . This leads to a bound \u p ' 1 ' l+l \ 12 < \u p ' l (l,.)\ 2 ( 1 + C— 

for some constants Cf, and C > 0. Next part (iv) of lemma 12.51 is used in 
order to estimate the 'correction terms' 5u p,k ' l+l with respct to the norm 
|.|o,i and |.|i,2- Note that we may choose C* dependent only on v and the 
dimension n large enough such that the bounds are independent of the time 
step I. 

3 Further remarks on uniqueness and regularity 

Uniqueness is a consequence of a global existence result for certain semilinear 
systems (cf. also proposition 8.6 below). 

Theorem 3.1. Let Q, = W 1 or £1 = T n . There is a unique globally bounded 
and Holder continuous classical solution of the Cauchy problem |7]) of the 
multivariate Burgers equation. 

Proof. Let v 1 and v 2 be two globally bounded and Holder continuous so- 
lutions of the Cauchy problem ([1]). The difference <5v := v 1 — v 2 with 
components 5vi := vj — vf satisfies 

{ dSvj v^n d 2 Svj , v^n „l*i iT" X„ - dv i — PI 
(112) 
6v = 0, 

Since vj and vf are given and Holder continuous, we have a semilinear 
system which is a special case of the linear system (2.8) described in ([I]). 
Since vj , vf are globally bounded and Holder continuous the assumptions 
(2.9)-(2.13) of proposition 2.3 in ([!]) are satisfied in case = R n and 
according to proposition (2.3) there is a unique global classical solution 
5v G C]' 2 ([0,T) x ]R n ,IR n ). Since Jv = is a solution of (fTT2l) . we have 
v 1 = v 2 . The result can be transferred to the case f2 = T n easily. □ 
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